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We investigate the impact of uniaxial strain on atomic shifts, dipolar interactions, polarization and 
electric permittivity in TiCh (rutile) by using two different implementations of density functional 
theory. It is shown that calculations using the Vienna ab inito simulation package (VASP) and the 
plane-wave self-consistent field method (PWscf ) yield qualitatively the same atomic relaxations and 
ferroelectric trends under strain. The phonon dispersion curves of unstrained and strained Ti02 
(rutile) obtained by employing the linear response method confirm previous calculations of the giant 
LO-TO splitting and the appearance of soft polar modes. A second order phase transition into a 
ferroelectric phase with polarization along (110) appears under expansive strain in (110) direction. 



I. INTRODUCTION 

TiC>2 is one of the key materials within oxidic sys- 
tems. Due to its large bio-compatibility and the well 
optimized and cheap processing, it is used in a huge va- 
riety of applications, e.g., see Tab. [TJ The first class of 
applications is based on its high photocatalytic power 
under UV-lightP^ 

Additionally, Ti02 is transparent for visible light and 
possesses a large permittivity and therefore high refrac- 
tion indices and is thus a very powerful material for opti- 
cal tasksP The large permittivity also gives rise to ap- 
plications of TiC>2 in microelectronics, for example as 
thin film capacitors. - Especially the rutile morphology 
of TiC>2 is important which is an incipient ferroelectric 
materialP Although, no ferroelectric transition has been 
observed experimentally until now, its dielectric constant 
in (001) direction is large and increases drastically with 
reduced temperature^ which indicates that the system 
is very close to a ferroelectric transition. The anomalous 
behavior of the dielectric constant can be linked to the 
softening of the optical A 2u mode, which is the only po- 
lar mode in (001) directionP Theoretical investigations 
have shown further softening of this mode under lattice 
expansion or strainP^ TiC>2 may be even interesting 
in the context of multifcrroics since magnetic proper- 
ties can be induced by doping with magnetic ions^ or 
when nanocomposite systems are formed. For example, 
magneto-electric coupling at Fe-TiC>2 interfaces has been 
observed experiment ally™ This would lead to new mul- 
tifunctional devices, see Tab. [Tj 

Most existing applications are based on rutile surfaces 
and interfaces, especially on the rutile (110) surface. This 
surface can be strained in (110) or (001) direction by the 
growth on a proper substrate. While a detailed discus- 
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sion of the influence of (001) strain and isotropic expan- 
sion exists, the influence of (110) strain is still not com- 
pletely understood. For example, the impact of (110) 
strain on the A2« mode is still under debate as this 
mode has been predicted to be insensitive with respect 
to such strain in Ref. [13l whereas we recently found a 
softening of the modeP^l Besides the destabilization of 
a "polar mode" in (001) direction corresponding to a 
strong ferroelectric tendency, other polar phonon modes 
in (110)/(1 10) di rection soften due to expansive (001) or 
(110) strainffSEH We found out that a stable ferroelectric 
state with a large polarization in (110) direction can be 
stabilized under expansive strain along (110) within our 
ab initio simulations. 12 Hence, stable polar states can be 
expected in both, (001) and (110) direction, which show 
different dependencies on the magnitude of the imposed 
strain and its direction. Also, a coupling between these 
soft modes exists.^ Because of this, the possibility of 
"strain engineering" the magnitude and direction of the 
polarization or at least of the dielectric and the piezo- 
electric responses, due to different relative magnitudes of 
imposed (001) and (110) strain, is very likely. 
Nevertheless, a detailed understanding of the strong po- 
lar tendency in (110) direction is still missing and our 
previous results are in part contradicting the theoreti- 
cal results of Ref. [T21 where the destabilization of the 
structure along an acoustic phonon mode, or a non-polar 
phonon mode has been predicted under (110) strain. 
The main goal of the present paper is to gain under- 
standing of the influence of (110) strain on the phonon 
spectrum and especially on the polar modes polarized in 
(110) direction. For this purpose, ab initio simulations 
have been performed with two different pseudopotential 
codes in order to cross-check our results. The atomic 
displacements and the ferroelectric properties have been 
monitored for different values of expansive (110) strain. 
Within both kind of simulations we find a ferroelectric 
instability of the system along the lowest polar mode in 
(110) direction. The phase transition to a ferroelectric 
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state with large polarization is of second order and its 
onset is below 0.5 % strain. Additionally, the lowest po- 
lar mode in (001) direction also softens, but to a lesser 
extent. 



II. COMPUTATIONAL DETAILS 

All results have been obtained by self-consistent den- 
sity functional theory calculations employing the plane 
wave pseudopotential codes VASP (Vienna Ab Initio 
Simulation Package pSl anc i the PWscf (Plane- Wave Self- 
Consistent Field) code which is part of the quantum 
espresso package.^ Within the VASP simulations the 
projector augmented wave potentials^ with the gener- 
alized gradient approximation (GGA) in the formulation 
of Perdew, Burke and ErnzerhoP2 for the exchange cor- 
relation potential have been used. For Ti (O) atoms 
4s3d (2s2p) electrons have been treated as valence in 
the calculations. Within the PWscf simulations ultra- 
soft pseudopotentials with the local density approach 
(LDA) and Ti (O) 4s3d3s3p (2s2p) valence electrons have 
been used. All ionic positions have been relaxed until 
forces converged to 0.001 eV/A. Born charges, the di- 
electric permittivity and phonon frequencies have been 
calculated using density functional perturbation theory 
(DFPT)P The "modern" theory of polarization^ has 
been used within our VASP simulation in order to deter- 
mine the spontaneous polarization. For the calculation 
of the phonon spectrum a Monkhorst-PaclP^ k-point grid 
containing 126 points was used leading to a convergency 
of at least 0.1 cm -1 for selected special points of the Bril- 
louin zone. The phonons were calculated using 4x4x4 
q-point grid. 



III. RESULTS 

The rutile morphology of Ti0 2 possesses P4 2 /mnm 
symmetry. Each Ti-atom is 6-fold coordinated by O- 
atoms with 4 short equatorial Ti-0 bonds and two longer 
apical bonds, see Fig. [T] Besides the lattice constant a 
and the tetragonal ratio c/a, an internal parameter u de- 
fines the relative length of both kinds of Ti-0 distances. 
The values, which have been obtained for these lattice 
constants agree with former results from literature, see 
Tab. |Tlj While the LDA lattice volume is slightly under- 
estimated in comparison to the experimental values, it is 
overestimated in GGA calculations. Although, this en- 
larged atomic volume may also lead to an overestimation 
of the ferroelectric character, the quantitative tr ends ca n 
be properly described by using GGA potentials! 9 * 11 * 27 * 

Previous work confirmed that the ferroelectric transi- 
tions of ruti le is linked to the softening of polar phonon 
modes ElSliSI Because of this, the phonon spectrum of the 
material will be discussed in the following. The irre- 
ducible representations of the optical vibrations of rutile 



at r are given by symmetry^ as 

r op t = Ai g + A 2g + A 2u + 2B\ U + Big + B 2g + E g + 3E U . 

Besides the lowest frequency polar mode, A 2u , the system 
possesses 6 modes of E u symmetry which are polarized in 
the (110) / (110) plane. Since all these modes are doubly 
degenerated, the eigenenergies E Ui i, E Ui2 and E Ui3 corre- 
spond to a basis of two eigenvectors each and thus atomic 
displacements in the whole plane are not distinguishable. 
One possible basis for the E Uji modes can be constructed 
from the eigenvectors given in Tab. |IV| combined with 
the same displacement pattern in (110) direction. 

The lowest mode, E Uj i, consists mainly of a rigid shift 
between Ti and O sublattices. Within the E Ui 2 mode the 
Ti-sublattices are shifted against each other hereby each 
Ti-atom is shifted into the same direction as its equatorial 
O neighbors. The highest frequency mode, E„ : 3, shows 
exactly the opposite trend with a relative shift between 
the Ti sublattices and the Ti-O e? neighbors, while the 
Ti atoms and their apical neighbors are shifted into the 
same direction. 

The double degeneracy of the E u modes at Y is lifted 
because the large dipolar interaction induces an en- 
ergy penalty for the longitudinal optical modes (LO) in 
comparison to the transversal optical modes (TO), see 

Tab. urn 

Besides the energetic LO-TO splitting, the dipolar inter- 
action induces a finite coupling between the E u modes 
and thus modified eigenvectors appear for the LO modes, 
see Tab. |IV| An alike mixing of the TO E u modes due to 
the dipolar interaction has been discussed in Ref. [5J 

Figure [2] shows the calculated phonon spectrum in 
parts of the first Brillouin zone. The eigenvalues ob- 
tained at r are compared to literature in Table |III| Our 
calculated phonon spectrum matches qualitatively pre- 
vious theoretical studies and experimental results. As 
in Ref. [13], we find an anomalously soft acoustic phonon 
mode in the Z-A-M region of the Brillouin zone. How- 
ever, it has to be noted that the absolute value of the 
corresponding frequency at this point is in the order of 
the accuracy of our calculation. Besides this, there exist 
large quantitative deviations. As the phonon frequen- 
cies are very sensitive with respect to small changes of 
the lattice constants, especially the lowest polar modesp^ 
these deviations can in part be explained by the slightly 
larger atomic distances obtained within our LDA calcu- 
lations in comparison to former results. Most notably, it 
has been showrP^ that the Bi g frequency increases with 
lattice expansion while the LO E„ modes remain nearly 
unchanged. Thus, our deviations of the mode frequen- 
cies are showing the expected trend for a small lattice 
expansion. 

Additionally, differences in the used basis sets, the im- 
posed sum rules and the accuracies of the simulation have 
drastic effects on the lowest modes. Detailed discussions 
for T1O2 and the similar incipient ferroelectric material 
SrTi0 3 can, for example, be found in Refs. H7J and [301 
The large quantitative deviations with respect to exper- 
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TABLE I. Devices based on Ti02 and future perspectives. 



Photo-catalytic properties Optical properties DMS J 'Electronic devices ME c and multifcrroic properties 



Solar ceP 
Self-cleaning materiaP 
Cancer therapy^ 
Air purificatiorP 
Sensors^ 



Anti-reflection coating Magnetooptic devices 

Pigments MRAM^] 

Interference filterP Spintransistor 

Optical waveg uideP . . . 

Sun protection Thin film capacitors^ 



Modulation of optical waves 
ME switching 
ME data storage 
Amplification 
Frequency conversion 



a Co doped Ti02 works as transparant diluted magnetic semiconductor (DMS) in Ref. 1141 

b A review on oxidic DMS materials can be found in Ref. 1161 

c See review \T7\ on magneto-electrical coupling (ME) and possible devices. 

d Magnetic Random Access Memory 
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FIG. 1. Atomic structure of Ti02 rutile. The approximate atomic size is fixed by covalent radii. Ti: light, cyan, O: dark, red. 
(a) Primitive unit cell, (b) V2 x y2 x 1 supercell. Ovals mark the different sublattices under (110) strain. Dotted line: class 
1; solid line: class 2, see text. Arrows represent the main atomic shifts under (110) strain for one subclass. 



TABLE II. Lattice parameters of bulk rutile. a: lattice con- 
stant; c/a: tetragonal ratio; u: internal parameter; ZX corre- 
sponds to the principal values of the Born tensor for Ti, along 
(001), along the apical Ti-O bond (||) and perpendicular to 
these directions. e^ i) an d is the permittivity along (001) 
respectively in the (010)/(100) plane. 

VASP PWscf RefP RefP Exp. 



a [A] 
c/a 
u 

"^(001) 

z \\ 
Zl 

oa 
e (001) 
oo 



4.664 4.581 4.536 4.572 4.59^; 

0.637 0.641 0.643 0.644 0.64^' 

0.305 0.305 0.304 0.304 0.30. 

8.15 7.69 7.54 7.77 

7.60 7.26 7.34 7.49 

5.36 5.26 5.43 5.43 

9.78 8.95 8.67 9.21 8.4. 

7.95 7.46 7.54 7.95 6.84 



b RefPH 



polar mode under lattice expansion or expansive 
(001) strain induces a ferroelectric transition. This effect 
can be explained by simple energy arguments. The rela- 
tive stability of the ferroelectric and paraelectric phases 
depends on the dipolar energy gained by the ferroelec- 
tric state and the energy penalty due to the short range 
repulsion which builds up for the relative atomic displace- 
ment. As the equatorial Ti-0 bonds contribute most to 
the short-range repulsion, a ferroelectric transition be- 
comes likely if this distance is sufficiently enlarged. Sim- 
ilarly, a unifor m l attice expansion causes a softening of 
the E u modespE^ although this effect is less pronounced 
than for the (001) direction. Additionally, one may ex- 
pect softening of the polar modes in (110) and (001) di- 
rection for expansive (110) strain in this picture. 



III.l. (110) strain 



iments can be partially traced back to the measurement 
at finite temperatures since it has been shown in Ref. EH 
that the lowest polar mode softens substantially with de- 
creasing temperature. 

As discussed in the introduction, the softening of the 



Besides the modification of the Ti-0 distances, strain 
in (110) direction reduces the symmetry of the system 
to the orthorhombic Cmmm class. Two inequivalent Ti 
and O sublattices form, see Fig. [TJb). Within the first 
subclass the short Ti-0 distances are conserved whereas 
the equatorial bond length within the second subclass is 
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TABLE III. Upper part: Frequencies of the optical phonon modes at F for undistorted rutile and the deviations between our 
PWscf simulations and experimental results. Lower part: Frequencies of the polar modes under 2 % expansive strain along 
(110). 





A 2u (TO) 


A 2u (LO) 


E u ,i (TO) E 


,,2 (TO) E„, 3 (TO) E 


u,v (LO) E u 


,v (LO) E u 


.3' (LO) 


RefP 


176.1 


769.3 


164.8 


391.3 


492.8 


351.5 


441.7 


808.4 


Exp. 


144* 


811 


183 


388 


500 


375 


429 


842 


RefP 


126 


749 


127 


375 


476 


342 


432 


782 


RefP 


47 




127 


357 


472 








PWscf 


38 


743 


90 


367 


475 


355 


429 


765 


Deviation[%] 


-73 


-8 


-50 


-5 


-5 


-4 





-6 



Bin (TO) Bu, (LO) bVi(TO) B 2ll , 2 (TO) B 2u , 3 (TO) B 2u , r (LO) B 2u , 2 ,(LO) B 2u , y (LO) 



i 175 i 72 157 362 487 358 439 753 

90 334 B 3 „,i(TO) B 3ll , 2 (TO) B 3u , 3 (TO) B 3M|1 >(LO) B 3 ^ 2 ,(LO) B 3m , 3 ,(LO) 

336 725 i227 327 454 325 385 684 




FIG. 2. Phonon spectrum of Ti0 2 rutile at the equilibrium lattice parameters obtained with PWscf (dotted lines) with 
highlighted polar modes at F (solid lines). Soft modes for 2 % expansive (110) strain are included (dashed lines). (1) Bi„ (LO) 
(2) B 3ll (TO) (3) Bi„ (TO) (4)-(5) acoustic (6) B 3ll (LO) at F. 



modified. The optical modes at T are now represented 

bySS 

r op t = 2A g + 2Bi g + B 2g + B 3g + 3Bi u + 3i?2u + 3i?3„ 

and the degeneracy of the E u modes is lifted to B 2U) j 
in (110) direction and B 3n ,i in the strained (110) direc- 
tion. An alike mode splitting appears for the acoustic 
modes and thus the destabilization of the atomic struc- 
ture alon g th e acoustic mode in (110) direction has been 
predicted 



Although, we find an alike softening of acoustic modes 
(branches (4)-(5) in Fig. [2]), this phonon branches show 
a large mode mixing and are both polar in (110) direction 
within the soft regions. Additionally, our previous inves- 
tigations with VASP yield a destabilization of the sys- 
tem towards a ferroelectric state under expansive strain 
caused by the softening of a polar modeP^ 
In the following, the atomic displacements and the re- 
sulting ferroelectric state will be discussed in more de- 
tail. Within these simulations no diagonalization of the 
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(110)-strain(%) (1 10)-strain (%) 

FIG. 3. (a) Atomic shift of the atomic subclasses in (110) direction under expansive (110) strain. Red squares: Tii; Blue 
diamods: Ti2; Cyan triangles: Oi; Black plus: 02,1; Green crosses: 02,2; The inset in (a) shows 0-1% strain, (b) Born 
charges of Ti within the ferroelectric state for different (110) strains. Tii: Red solid curve, Ti2: Black dotted curve; Diamonds 
Z(ooi)' Triangles: Z*|j Squares: Zj_. (c) Spontaneous polarization for different values of strain, (d) Evolution of the electronic 
permittivity within the ferroelectric state with strain. Black diamonds: (001) direction; Red squares: (110) direction; Blue 
triangles: (110) direction. 



dynamic matrix has been performed. Instead the para- 
electric symmetry of the system has been broken by a 
shift of one O-atom within a y/2 x x 1 supercell, see 
Fig. [l] The atomic positions have been optimized within 
this Pmm2 symmetry class and thus, the O2 atoms are 
split into two inequivalent subclasses. 
The amplitude of the obtained polar shift increases con- 
tinuously with imposed expansive strain, see Fig. |3][a). 
Already at 0.5 % strain, the amplitude of the displace- 
ment along the polar mode is exceeding the arbitrary 
displacements below the given accuracy of the atomic re- 
laxation. Between 2 % and 6 % strain, the displacements 
of Tii and 02,2 atoms increase approximately linearly, 
while the shift of the Ti 2 and Oi atom is nearly constant 
within this range of strain and the increase in the 02,1 
displacement is only minor. For larger values of strain 
the Tii amplitudes start to saturate, while the Oi and 
Ti 2 displacements even decrease. This evolution of the 
atomic displacements with imposed strain is a fingerprint 
of a second order phase transition. Although, the ex- 
act transition strain cannot be determined, the transition 
seems to start for infinitesimal expanding strain. Most 
notably, the spontaneous polarization shows the charac- 
teristic profile of such a phase transition with a continu- 
ous onset and a saturation for large values of strain, see 
Fig.§c). 



The saturation of the spontaneous polarization can 
again be explained by simple energy arguments. While 
the reduction of the short range repulsion further in- 
creases with large values of strain, the dipolar energy is 
reduced in this case. As Born effective charges are pro- 
portional to the dipolar interaction, this argument can be 
confirmed by the evolution of the Born charges in (110) 
direction, Z^ 110 ^, under strain, see Fig. |3j[b). The ten- 
sor of the Born charges is diagonal within a (110), (010), 
(001) reference system, with the principal values Z^ 001 ^ 
aligned along (001) direction, Z|j aligned along the apical 
bond and Z* ± perpendicular to both other directions. For 
a detailed analysis, see Refs. ITT1 and 1521 

For Tii atoms, the apical Tii-02 bonds are aligned 
along (110) and thus the large dipolar interaction leads 
to a large shift of the corresponding Ti-atoms, cf. Z|j. 
Instead, the apical Ti 2 -Oi bonds are aligned along (110) 
and thus the dipolar interaction in (110) is much smaller, 
cf. Zj_. For both Ti classes, the Born charge along (110) 
decreases considerably with increasing strain as does the 
dipolar interaction. Especially, for the Ti 2 atoms the 
Born charge is reduced to the nominal charge of Ti 4+ at 
about 6 % strain and thus vanishing dipolar interaction 
and atomic shifts appear for increasing strain. 

In order to characterize further the properties of the 
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TABLE IV. Displacement pattern of the mode eigenvectors 
(arbitrary units). For each atomic class, see Fig. [Tj the am- 
plitudes of the atomic shift in the given directions are listed. 
Relative amplitudes of the B u ,; modes are given for 2 % ex- 
pansive strain. 







Tii 


Ti 2 


Oi 


o 2 


E Ul i/B2u,]| 




0.1 


0.3 


-0.3 


-0.3 


E u ,2/B2ii,2 




0.3 


-0.2 


0.2 


-0.4 


E u ,3/B2ii,3 




0.2 


-0.3 


-0.2 


0.4 


E„ l i'/B 2 „ l i' 


(LO) 


0.2 0.0/-0.1 0.2 


-0.4 


E u .2'/B2ii,2' 


(LO) 


0.3 


-0.5 


-0.1 


0.3 


E u ,3'/B2u,3' 


(LO) 0.2 


0.2 


-0.5 


-0.1 


B 3u ,i(T0f 




0.3 


0.1 


-0.2 


-0.4 


B 3 „,2(TO) 




0.1 


-0.3 


0.4 


-0.2 


B 3 „, 3 (TO) 




0.2 


-0.3 


-0.2 


0.4 


B3„,i'LO 




0.2 


-0.3 


0.4 


-0.2 


B3„,2'LO 




0.4 


-0.2 


-0.4 


0.1 


B3„,3'LO 




0.1 


0.2 


-0.1 


-0.5 



a Displacement in (110) direction 
b Displacement in (110) direction 



Tii Ti 2 Oi 2 



A 2u (LO)/(TO 
Bi„,i(TO) 
Bi„, 2 (TO) 
Bi«, 3 (TO) 



0.3 0.3 0.5 0.5 

0.6 -0.5 -0.3 

0.7 -0.3 -0.3 -0.4 

0.1 0.5 -0.5 



Bi„,i(LO) 
Bi„, 2 (LO) 
Bi„, 3 (LO) 



0.7 -0.7 0.1 -0.1 
0.3 0.3 -0.4 -0.5 
-0.1 0.0 0.5 -0.5 



Displacement in (001) direction 



ferroelectric state, its electronic permittivity, Coo, for dif- 
ferent values of strain is shown, see Fig.[3]jb). Although, a 
large permittivity is obtained over the whole investigated 
interval starting form % to 11 % expansive strain, all 
principle values of the permittivity tensor decrease more 
than linearly with increasing strain. Especially, in 
(110) is considerably reduced. 

Additionally, the influence of the strain induced ferro- 
electric state on the electronic structure has been inves- 
tigated. For this purpose, the density of states has been 
analyzed for the ferroelectric state at 5 % strain. As only 
small deviations in comparison to the undistorted mate- 
rial appear, similar electronic properties have to be ex- 
pected. Only, the obtained band gap is slightly enlarged 
to 1.82 eV in comparison to the band gap of 1.64 eV for 
the undistorted bulk. 

If the polar shift is prevented by the imposed Cmmm 
symmetry, or due to the large short range repulsion which 
appears under compressive strain, only the oxygen atoms 
relax and thus the internal parameter u is modified. This 
relaxation corresponds to a change of the relative Ti-O e(? - 
Ti-O ap distances due to the larger compressibility of the 
apical bond in contrast to the equatorial bondP Hence, 



the apical bond length is further enlarged for expansive 
strain, while the strain-imposed increase of the equato- 
rial bond length is reduced (vice versa for compressive 
strain). Qualitatively the same relaxations are obtained 
within VASP and PWscf simulations. The main relax- 
ation can be linked to a Ai g mode of the undistorted 
system for which we yield a frequency of 588 cm _1 at 
r. Since the strain induced modification of the Ti-0 dis- 
tances is not uniform in (110) direction and (110), the 
relaxation of the internal parameter differs in both direc- 
tions and thus a modification of the relaxation amplitude 
for the O-subclasses is superimposed. 

In Ref. [12] we have shown that the calculated ferro- 
electric displacement pattern can be approximated by a 
simplified mixture of the bulk E u modes with unique am- 
plitudes. In agreement with this finding, the diagonal- 
ization of the dynamic matrix at 2 % expansive strain, 
based on PWscf, yields an imaginary phonon frequency 
for the lowest mode. The eigenvector, see Tab. IV 
of this mode is qualitatively matching our earlier results 
obtained with VASP, as the main shift consists of a rela- 
tive shift between the Ti and O atoms with superimposed 
relative shifts between the Ti and O atoms corresponding 
to different subclasses. Note, that the PWscf results have 
been obtained for the primitive cell wit Cmmm symme- 
try in order to reduce the computational time and thus 
the O2 atoms are all symmetrically equivalent. There- 
fore, the appearance of the discussed ferroelectric state 
with (110) polarization can be qualitatively reproduced 
for different simulation details. 

Besides the large softening of the lowest ^ u mode, all 
other B 3ll modes are shifted to lower energies at T due to 
the discussed lifting of the degeneracy. In contrast, the 
B2 U modes are shifted to higher frequencies. The only 
exceptions are the ^>2u,2 and ^2u,y modes which are also 
lowered in comparison to the undistorted bulk. 

The strain imposed modification of the eigenenergies is 
accompanied by a modification of the mode eigenvectors 
for the B% u modes, while the B^ u eigenvectors are iden- 
tically to the corresponding E u modes of the undistorted 
system within the given accuracy. 

The relative amplitudes within the polar modes can 
again be qualitatively understood by the counterplay of 
dipolar interaction and short range repulsion. As pre- 
viously discussed for the (110) direction, the Ti atoms 
belonging to class 1 contribute more to the dipolar in- 
teraction than the Ti2 atoms. The opposite trend holds 
for the shift in (110) direction, where a larger contribu- 
tion to the dipolar interaction is obtained for the second 
subclass. Additionally, the equatorial Ti2-C>2 distance 
is enlarged due to the imposed strain and thus the short 
range repulsion is reduced for a relative shift between the 
Ti 2 and O2 atoms. Both effects lead to a larger Ti 2 -02 
shift in comparison to the Tii-Oi shift for the B2 U modes. 
Although, a similar short-range repulsion is obtained for 
the B3„ modes, the different dipolar interaction leads to 
modified relative amplitudes. 

Besides, the discussed softening of the B3„ modes, our 
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previous investigation yields a softening of the A 2li mode 
under expansive (110) strainj^if the polar displacements 
in (110) direction are prevented by the imposed symme- 
try In contrast to this, the mode has been predicted to 
be quiet insensitive towards (110) strain!^ 
For the strained system the single polar mode in (001) di- 
rection, the A 2u mode, is replaced by 3 polar B lu modes. 
The corresponding mode eigenvectors are listed in the 
lower part of Tab. [TV] While the A 2u LO and TO modes 
are characterized by a shift of the Ti-0 sublattices with 
a unique amplitude, relative shifts between the Ti- and 
O subclasses are superimposed for the Bi u (TO) modes. 
In Ref. [T^lwe assumed unique amplitudes for the atomic 
displacements and thus the atomic shift along the real 
mode pattern is most likely even more favorable. 
The diagonalization of the dynamical matrix for 2 % 
strain yields imaginary frequencies for the (LO) and (TO) 
Bi„.i modes. Therefore, our results obtained with the 
PWscf code confirm our previous findings that the ferro- 
electric state polarized in (001) direction is energetically 
more favorable than the paraelectric state although less 
favorable than the ferroelectric state polarized in (110) 
directionPSl 

It is most likely that a ferroelectric transitions along (001) 
takes place along the Bi„ (TO) eigenvector. For this 
eigenvector, the relative shift between the Ti 2 -02 equa- 
torial bonds is enlarged in comparison to the shift of the 
Tii-Oi neighbors, due to the different short range repul- 
sion. If the evolution of the dipolar interaction is again 
monitored by the Born charges, an even larger modula- 
tion of the atomic shifts has to be expected for increas- 
ing strain, as Z^ 0Q1 ^ is decreasing for Tii with increasing 
strain, while the dipolar interaction is nearly constant for 
the Ti 2 atoms. 

As in SrTiOjP^ a large coupling between the soft modes 
in (110) and (001) directions exists and the atomic dis- 
placement along the B lM modes is disabled within the 
polar state in (110) direction. This means, the phonon 
frequencies for an polar displacement along (001) are no 
longer imaginary if the atoms are displaced along the B3 M 
eigenvector. 



IV. CONCLUSIONS AND OUTLOOK 

In this paper we have discussed the softening of po- 
lar modes under (110) strain in rutilc, finding a second 
order phase transition (polarization versus (110) strain) 
for a polar state in (110) direction. The calculations were 
done using the plane-wave method as implemented in the 
VASP and PWscf codes. The two methods allowed for 
a cross-checking of the results and lead to qualitative 
agreement regarding the atomic relaxations and atomic 
mode patterns under (110) strain. The decomposition 
of the displacement pattern "under strain" into ideal E u 
modes is qualitatively matching the eigenvectors of the 
lowest i?3 M mode and a detailed analysis of the polar 
modes and associated eigenvectors has been undertaken. 
The investigation of displacements, dipolar interaction, 
dielectric permittivity and electronic structure underlines 
the trend for a second order ferroelectric transition un- 
der (110) strain in rutile Ti0 2 . The calculations have 
been complemented by evaluating the phonon dispersion 
curves of unstrained rutile and the appearance of imagi- 
nary branches at 2 % strain, which correspond to polar 
displacements in (001) and (110) direction. The coupling 
of these polar B\ u and B^ u modes as well as their dif- 
ferent dependencies on (001) and (110) strains leads to 
the possibility of strain engineering. Thus, a phase dia- 
gram for the relative stability of both kinds of modes at 
different combinations of strains has to be calculated in 
future. In addition, it will be necessary to calculate the 
phonon spectra of both ferroelectric phases and to char- 
acterize their interrelation. This should be supplemented 
by evaluating the electron-phonon coupling constant and 
its impact on the structural transitions, which may oc- 
cur in strained rutile. Finally the relation to experimen- 
tal setups has to be achieved. Therefore, calculations on 
strained surfaces will be performed in order to investigate 
the influence of the surface and surface induced atomic 
relaxations on the polar modes in thin films of Ti0 2 . 
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